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Abstract 

The form factors for the Ki^ semileptonic decay are computed to order O(p^) in 
generahzed chiral perturbation theory. The main difference with the standard O(p^) 
expressions consists in contributions quadratic in quark masses, which are described by a 
single divergence- free low-energy constant, A^. A new simultaneous analysis is presented 
for the CKM matrix element Vus, the ratio Fk/F-,^, K^^ decay rates and the scalar form 
factor slope Aq. This framework easily accommodates the precise value for Vud deduced 
from superallowed nuclear /3-decays. 
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1 Introduction 



Together with low-energy tt-tt scattering |jT|, ^ , the semi-leptonic i^^s and i^e4 @, H] decays 
have been among the first apphcations of standard chiral perturbation theory (S^PT) [|], 0, |^ 
to be studied at the one-loop level. At this order, the mesonic form factors that describe these 
decays do not contain the poorly known low energy constants Lg and Lq, and consequently 
they may be expected to be less sensitive to the size of the chiral condensate (qq) than, e.g., 
the TT — TT s-wave 0, 3. In particular, the K^^ form factors at O(p^) merely involve, besides 
the masses and decay constants of the pseudoscalar states, the low energy constant Lg, whose 
value can be obtained from the experimentally measured charge radius of the pion [^. This 
fortunate circumstance has been used in the past to extract the CKM matrix element V^s from 
the Ke3 decay rates and presently this extraction is still considered by the Particle Data 
Group (PDG) compilation |]12| to remain the most accurate and least model dependent. Yet 
this determination of Vus relies on a mo del- dependent estimate of O^m'^^^j.^J contributions to 
the form factor /+(0) (the notation will be given below), which, in S^PT, arise as contributions 
of order O(p^). With the new careful and accurate determinations of Vud from superallowed 
nuclear /^-decays fl^, [1^, |1^ , the size of these corrections is required to be comparable to the 
(parameter-free) genuine 0(p^) contribution in order to preserve the unitarity of the CKM 
matrix. It is then legitimate to ask whether similarly important Olm^^^j.^) contributions would 
not affect the parameter-free S^PT prediction for the slope Aq of the K^^ scalar form factor 
at order 0{p^). In principle, these questions can be answered by performing the full two-loop 
SxPT calculation of the Ki^ form factors, provided the several new 0{p^) counterterms that 
will contribute can be measured independently or estimated in a reliable way 0. The present 
status of this enterprise is limited to the evaluation of the so-called chiral double-logs |jl9 



which, although only part of the full two-loop corrections, do not seem to point towards huge 
effects coming from the chiral loops themselves (see in particular the numerical results in Table 
1 of Ref. [0 and the comments preceding it). 

The purpose of the present work is to address some of these issues from the point of 
view of generalized chiral perturbation theory (G^PT) [§, ^ The expressions for the 
form factors of semileptonic decay of kaons to O(p^) in G^PT have never been published (a 
discussion of the Ki4^ form factors at order O(p^) in G^PT can be found in Ref. |2^), although 

^The structure of the order 0{p^) effective lagrangian of SxPT has been discussed in Refs. 
^Unfortunately, the 0{p^) calculation of Ref. ||l^ only considers a very specific combination of the slopes of 
the /J^'^ and of the mesonic electromagnetic form factors. 
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they have existed for some time and have been partially reported on various occasions. Keeping 
in mind the well-known and important exception of the phase of the K^^ decay amplitude, the 
and Kg^ form factors are indeed found to be to a large extent independent of the size of 
the condensate (qq). The main consequence of the modified chiral counting [rriq ~ C^(p)) is 
that in G^PT the form factors receive a contribution quadratic in quark masses already at 
order 0(p'^), in addition to the standard 0{p'^) expressions. Furthermore, these new C^(?^quark) 
contributions, which within S^PT would count as 0{p^), are all related. They all stem from 
a single term of the ^(^2,2) component of the effective Lagrangian Ces and they are described 
by a single divergence-free low-energy constant A3 (see Ref. ^ and Appendix A for notation). 
(This statement is exact in the case of Ki^ form factors /+ and /_ , whereas in the case of Kf,^ 
it remains true for the dominant 0{m1) terms.) It thus appears that the O(p^) G^PT offers a 
predictive description of the terms quadratic in quark masses which is non-trivial compared to 
the one-loop S^PT (no quadratic terms) and yet much simpler than the standard 0{p^) order, 
in which other unknown constants should contribute in addition to the 0{p^) G^PT terms 
driven by the constant ^3. This framework, which is as systematic as the standard expansion, 
suggests a new simultaneous analysis of Vus, Ki^ decay rates, Fk/F^,, and of the scalar slope Aq 
which may be of interest in connection with the constraint of unitarity of the CKM matrix and 
with the forthcoming new data. A closely related application, namely the determination 
of the C(p^) constant L3 from the K'^^ decay rate, will be presented elsewhere ||2^, together 
with a full discussion of the Ki4 form factors at order 0{p'^) in G^PT. 

The present stage of our analysis involves one additional limitation, to the extent that no 
electromagnetic corrections are included, unless explicitly stated. For this reason we postpone 
a detailed analysis of the isospin asymmetry in the K^^ and -K"^ decay rates that is due to the 
mass difference 171^ — rriu. We just check that this asymmetry is consistent with the G^PT 
treatment of 7r° — 77 mixing within errors. 

The paper is organized as follows: Section 2 provides the necessary expressions of the Ki^ 
form factors and of the 7r° — 77 mixing angles in G^PT. Implications for the determinations of 
Vus and of the ratio Fk/F^, of pseudoscalar decay constants are discussed in Section 3. The 
slopes of the K^^i form factors are considered in Section 4. Concluding remarks are presented 
in Section 5. Details on the 0(p^) structure of the G^PT effective Lagrangian and on its 
renormalization are presented in Appendix A. Useful expressions for the pseudoscalar masses 
and decay constants have been gathered in Appendix B. 
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2 Kis form factors in G^PT to 0{p^) 



We consider the two semileptonic decay channels 



K\p) 



TT {p)t{pi)yi{Pv 



a]- 



(2.1) 



The symbol £ stands for /i or e. As stated before, we do not consider electromagnetic corrections. 
The processes are then described by four form factors, and f^°'^ (t), which 

depend on t = (p' — pY = {pi +Pu)'^, the square of the four momentum transfer to the leptons, 
and which are defined in terms of the hadronic matrix elements of the charged strangeness 
changing QCD vector current as follows: 

1 



<n\p')\{s^,um\K^{p)> 

<n-ip')\is^,u)mK%p)> 



V2 



(2.2) 



2.1 GxPT expressions of the Kg^ form factors 

The one- loop G^PT expressions (using the notation of reference [Q) for the form factors are 
summarized below. We start with the two form factors f^°^ (t) and f^^^°{t), which are in 
practice sufficient for the description of the electron decay modes K^^ and K'^-^, and keeping 
isospin breaking contributions due to the quark mass difference — m„. For the K^^ channel, 
which is somewhat simpler, since — rj mixing only enters the loop contributions, we obtain 



r]K- 



it) 



+ V3i 
■1 



+ 



m 4 + -m A3 



(2.3) 



with (the definitions of the loop functions Mpg (t), Kpqit) and Lpq (t) that we use below can 
be found in Ref. 0) 



HpQ{t) 



1 



tM'pQ{t) - Lpgit) + -Lit 



(2.4) 



Here e denotes the leading order n'^ — rj mixing angle. 



where 



ITLs/m 



r2 = 2 



M 



^GMO 



K_ 

-2 



- 1 



(2.6) 



"iMl - AMI + Ml . 



We work only at first order in the quark mass difference {rrifi — rriu) , i.e. we consider only terms 
that are at most of order 0{1/R) where 



R 



— m 



(2.7) 



The lowest order S^PT value for the tt^ — rj mixing angle, £st= 4^, is recovered by dropping, 



4_R' 



in Eq. ( |2.5| ), the last two terms, which are counted as order O(p^) in S^PT. The and 
terms in Eq. (|2.3|) are 0{p^) contributions in G^PT (see the detailed formulas for the effective 
Lagrangian in Appendix 0), but are absent at this order in S^PT. 



For the KX decay mode, we find 



pK°tt 



V3, 



(2.8) 



v3 



Li? ^ 73. 



where the ir^ — r] mixing angles at order (9(p^), £1, £2, are defined in section |2^ below. Although 
the constant C2 corresponds to a counterterm of C(^2,2) (see Eq. ( [A.10| ) in Appendix 0) that 
violates the Zweig rule, it is not expected to be suppressed, since this violation occurs in the 
0~ channel. 

At zero momentum transfer, Eq. ( p.3| ) gives 



/f^"(0) 



1 



1 



2567r2F2 

r2 I ;\ #2 



2 (M2+ + M|o) ho' 



-Ml 



Ml,' 



+ [Mi, + Ml^)ho{j^) + ?>{Ml^ + M'^)ho{ 



K+ 



M2 



+ 2v^£ 



/\#2 /M'2 • 

{Ml + Ml) ho{^) - {Ml + Ml) ho{^) 



•Mi' 



M2 

rj 



+ 



8^2'^. 



r 



^We have however kept the pion mass difference, which is mainly an electromagnetic effect 
loop contributions. 



(2.9) 



251, in the 



4 



with 



ho{x) 



1 + 



2x 



\nx . 



(2.10) 



1 - 

To recover the expression at order 0{p'^) in S^PT [|n|, ^ for /|^°'^ (0), one replaces e by the 
leading order S^PT value £st, and one simply drops the last line in Eq. (|2.9|). Notice that 
this last contribution is the only correction of order 0{p'^) that does not vanish in the large- A^c 
limit of QCD. The fact that the O(p^) corrections to f^°'" (0) in S^PT vanish altogether in the 
large- Ai'c limit might provide a natural explanation why contributions of the 0{p^) counterterms 
could be comparatively sizeable. 

The expressions of the form factors /5"^ (t) and f^^^\t) are rather cumbersome if 
0(md — rriu) effects are included. Since we do not need the latter in this case, we give only the 
common expression of these two form factors in the isospin limit, which is 



= -m^(r-l) - -m^^^(r- l)(r + 3) - m^^^(r - l)(r + 2) 
+ im2(Ai + 5i)(r2-l) + ^m\A2-2B2)ir-l) 



+ m^D^{r~l){r + 2) + - 2/ii^ - 3^^] 

1 



(2.11) 



4Fi 
1 

4Ff 



[5t-2M^-2Mf^ + 16m'{Ao + 2 Z^){r + 1)] K^Kit) 
[3t - 2Ml - 2Ml + 4m2(Ao + 2Z^){r - l)(r + 3)] K^^(t) 



Mi 



M'k. it) + -LI 



+ 



For the ease of comparison, we may rewrite this G^PT expression for the form factor /_ {t) in 
terms of the corresponding S^PT expression: 



- 1 



2 LI {Ml -Ml) 



' F^ 
{2Ml + 2M^-3t) 



4F2 

{2Ml + 2M^-5t) 
3 (Mi - M2) 



K^K{t) 



2F2 



(2.12) 



j:K7T]std 

_\(Fk_ 
4 



K 



+ 2- 



F^ 



- 1 



5 



^m^ A3 (r- l)(r + 3) 

Agmo (r + 3) M,^ (1 + r) e(r) 



where 



Am'^ (Aq + 2 Z^) r2-r, 
- ^ = 2;l3l(l + 2C) 

C = ^oMo, (2.14) 

and the C(p^) expression 

-(r + 1) (Ai + 5i) - {A2 - 2B2) - (r + 3) A3 

-2{r + 2)m^D^ (2.15) 

+2m^^^ + 2(r + 2)m2^f 
1 

- (5/i^ - 2/ii^ - 3yu^) 



2(r- 

has been used. It is remarkable that, in Eq. ( 2.13| ), all ^^(2,2) constants - except for A3 - have 
been absorbed into the renormahzation of the decay constants. Notice also that, in contrast to 
/+(t), f-{t) starts only at next-to-leading order in the chiral expansion, and that furthermore 
this contribution is not suppressed at t = ioi Nc ^ 00. 



2.2 Mixing angles 

The expressions of the Kg^ form factors given in the preceding subsection involve the vr" — ry 
mixing angles e, 61,62- In defining them we ignore isospin breaking through electromagnetic 
effects, so that the only source of isospin violation is the quark mass difference — If 
7^ rnu, the isosinglet and isovector axial currents A^ and have nonvanishing off-diagonal 
matrix elements between the vacuum and one-meson states. The two n^r] mixing angles 61 
and 62 are introduced such as to define combinations of the axial currents having vanishing 
off-diagonal matrix elements: 

coseiAj(O) -sineiAj(O) |7r°(p)) = 

(2.16) 

{^\ cos 62 Al{Q) + sin 62 Al{Q)\r]{p)) = 0. 
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Both £i and £2 are of order 0{md — itlu). Note that we were not forced to define mixing 
angles in terms of matrix elements of the axial currents: we could have chosen to use the 
pseudoscalar densities, or any other operators with the appropriate quantum numbers; however, 
the corresponding expressions for the mixing angles would in general differ from those given 
below. 

Keeping only contributions which are at most linear in the quark mass difference md — rriu, 
the off-diagonal matrix elements of the flavor neutral axial currents read 



(Q|4(0)|7r°(p)) = ip.F^s^ 



(2.17) 



The decay constants and define the diagonal matrix elements of the same currents, 

(l]|4(0)|7r°(p)) = tp,F^ 



(2.18) 



Up to corrections of order 0((m„ — m^)'^) that we neglect, F^r can be identified with the charged 
pion decay constant, F-^ ~ 92.4 MeV. 

At order C(p^), the two mixing angles coincide and are equal to £. At order C(p^), they 
both receive C(mquark) corrections, and are no longer equal. Explicitly, 



£1 + £2 o 
^ = £ + 



v^i?(M2 - Af2) 



TT 



- 1 



^GMO , 1 / 
W 2 — r 

2 3^ 'Vr + 1 



f2 



and 



£1 - £2 



V^R{Ml-Ml) \F^ 



2, ,fr-2 



(2.19) 



(2.20) 



where we have defined 



X,{r) = (r - 1)^ [(3pi + 4p3) (r - 1) - P2 (r + 1)] 



(2.21) 



in terms of the >C(o,3) quantities pi,2,3- 
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3 Decay rates, Vus and Fk/Fj^ revisited 



3.1 The 



TT e^v rate 



For the i^es decay we only need to consider the form factor f+{t), and for K'^ Ti~e^v Eq. ( |2.3|) 
is all we need. Ignoring the tiny term in Eq. ( p.3|) proportional to we may write the decay 
rate in units of 10~^^ MeV as 



T[K^ ir-e+u] = V^, 105.056 

+ 203.031 (r-1) 
+ 98.267 (r - 1)^ 



2 f 2 X 2 



+ 



(3.1) 



Strictly speaking, the last contribution on the right-hand side of this expression, although 
coming from the square of the O(p^) expression of /^"'^ (t), represents an order O(p^) effect 
in the chiral expansion of the decay rate. However, as can be checked e.g. on Figs. 1 and 2 
below, it does not affect our analysis in the range of values considered for the Cabibbo angle. 

We express in terms of Vus as follows: first, we use Eq. ( p.l5| ), but only at lowest 
order, 



K 



- 1 



r - 1 V F2 

and then use the formula for the ratio of the branching rates 

TjK ^ H ^ Fl Mk+ {1 - {MjMK+ff 

r[7r ^ H {1 - {MjM^^fY 



(3.2) 



{1 + 5k- 5^) 



(3.3) 



where the radiative corrections 5k = 0.0020 ± 0.0002, = 0.0017 ± 0.0002 nearly cancel [gg 

Next, we use the experimental values |]T2| for the branching rates to fix ||TT|] the combina- 
tion of constants 

Fk Vus 



0.2758 ±0.0005, 

Ftt Vud 

together with the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 



(3.4) 



(3.5) 



We ignore \Vub\i which has been obtained [^, ^ from a model-dependent analysis of data 
from B semileptonic decays to be (3.25 ± 0.6) xlO"'^; however, even if the central value of Vub 
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turned out to be three times larger, this would not affect our results. After expanding Eq. ( |3.1| ) 
in Vus around Vus = 0.220, this implies 



5.390 + 31.10 0.220) 

+ (r - if m^As [5.059 + 37.58 (Ks - 0.220)] 

+ (r - l)^m'^Al [1.189 + 10.81 (Ks - 0.220)], 



(3.6) 



to be compared with the experimental rate of 4.937 ±0.052. Consequently, from the knowledge 
of Vus one may extract a value for the ^^(2,2) quantity m^(r — ifA^ = {uis — rfiYA^. For 
example, in Figure |l] we show the band in the {Vus , ^3) plane indicated by experiment, together 
with lines of constant Vus corresponding to the Particle Data Group (PDG) |]T2| value Vus = 
0.2196 ± 0.0023. Accepting this constraint on Ks would imply m'^{r - ifA^ = -0.089 ± 0.024. 
A naive dimensional analysis (NDA) [^ would give 



l)'m'A3\ 



NDA 



^s(Ah) 



0.04 



(3.7) 



where we have taken itLsIAh) ~ 200 MeV, and ^ 1 GeV is a typical hadronic mass scale. 



0.23 



0.225 



Vus 0.22 



0.215 




-0 . 15 



. 05 



(r-l)^ A3 



Figure 1: The dependence of Vus on A^. Horizontal lines indicate the range of values for Vu. 
given by PDG [|12I. 



3.2 Determination of Vud from nuclear beta decay 

In its 1998 update, PDG [|l^] recommends for Vus only the value determined from Kg^ decay, 
arguing that the value obtained from hyperon decays (see [^ for an early attempt along these 
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lines) suffers from theoretical uncertainties due to first-order SU(3) symmetry-breaking effects 
in the axial-vector couplings. However, within the context of the present G^PT analysis, we 
may not, without independent knowledge of ^3, use the values for \4s thereby extracted from 
Ke^ decay. Moreover, Fk is determined from the K^2 decay together with the knowledge of 
Vusi and we note that Fx/ F^^ appears implicitly in our theoretical expressions for the K^^ form 
factors. The origin of this is, of course, our re-expression Eq. of the >C(o,3) constant ^ in 
terms of Fk/Ft,. 

iFiom. the unitarity condition for the CKM matrix, it follows that \Vus\ may be fixed from 
the knowledge of the up-down quark- mixing matrix element of the CKM matrix, \Vud\-, alone. 
The value of Vud can be determined from several independent sources: nuclear superallowed 
Fermi beta decays, free neutron decay, and pion beta decay. 

Currently, superallowed Fermi 0+ 0+ nuclear beta decays [0, |l5l, together with 
the muon lifetime, provide the most accurate value. 



ud 



0.9740 ± 0.0005. 



(3.8) 



The precision is limited not by experimental error but by the estimated uncertainty in theoreti- 
cal corrections In Figure^ we show the band in the Vud , ^3 plane indicated by experiment, 
together with lines of constant Vud corresponding to these values. 



0. 977 



. 976 



Vud 0.975 



0. 974 




Figure 2: The dependence of Vud on ^3. Horizontal lines indicate the range of values for Vud 
determined from superallowed nuclear beta decay. 



The determinations of Vud from free neutron decay data are approximately a factor of 



four poorer in precision (see the discussion in |jT5[ and references therein), 

Vud = 0.9759 ± 0.0021, 



(3.9) 
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due to the difficulty in separating the vector from the axial piece, but planned experiments 
aiming at an accurate measurement of the electron emission asymmetry |^ could change this 



situation qualitatively, since the error in this case is primarily of experimental origin. 

Finally, the theoretical corrections in the nuclear Fermi transitions are absent in the case 



of the pion beta decay. The present status of this type of experiments results in the value [15 



= 0.9670 ±0.0161. (3.10) 



Here also, the situation might improve in the future ^ |0|. 
Combining the above results, one obtains 

Kd = 0.9741 ±0.0005. (3.11) 

Accepting this value for Vud would imply m?{r — Vj^A^ = — 0.124±0.022, somewhat larger than 
the NDA estimate Eq. (p.7|), but still acceptable. 

Finally, we note that recent (model-dependent) analyses of hyperon semileptonic decays 
give H 

Kd = 0.9750 ± 0.0004, (3.12) 



and 33 



Kd = 0.9743 ± 0.0009. (3.13) 



Taking the value Eq. (|3.11|) for Vud {i-e., excluding results from hyperon decays), the 



unitarity relation Eq. ( |3.5| ) gives 

Vus = 0.2261 ± 0.0023. (3.14) 

Incorporating the above values for the CKM matrix elements into Eq. (|3.4|) then implies 

Fk 

— - = 1.189 ± 0.012 [Vud from nuclear beta decay], (3.15) 

which may be compared with the corresponding result using the PDG values for Vus, 

Fk 

= 1.226 ± 0.014 [Vus from PDG]. (3.16) 

Using Eq. (p.4|), which relates Fk/F^, to Vus and Vud, Eq- ( |3.5| ) which relates Vus and Vud, and 
Eq. (|3.6| ) which relates Vus and (r — 

1)'^3, we may directly relate Fk/F^^ and m^(r — 

see Figure |. 
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-0.2 -0.15 -0.1 -0.05 

(r - 1)^ A3 

Figure 3: The dependence of FkIF-k on Aj,. 



3.3 mu 7^ rad effects on the n^e^v rate 

We may treat the decay of the similarly; the result is most conveniently expressed in terms 
of the ratio 



T[K^ 7r°e+z/] 



{R 



0.49615 - 0.00240 (r - 1 



2 /m^As m^^^ 



+ 0.00234 (r - 1)^ ( + — ^ 



V 2 



(3.17) 



where R'^^ is the ratio 



R^' ^ T^To-B ■ (3-18) 



/r"°(o) 



Note that the last two terms in the brackets in Eq. (|3.17|) arise from the small differences in 



phase space which are due to mass differences between and K~^, and between vr"*" and 7r°. 
Proceeding as we did for the rate above, i.e., using Eqs. ( p.2| ) and (p. 41), we obtain 

= (/?+0)2 {0.4961 + 0.0040 (Ks - 0.22) 

-(r - if A3 [0.0011 + 0.0041 (Ks - 0.22)] 

+0.0005 (r - l)^mM^} . (3.19) 

The correction terms in the curly brackets are completely negligible for A3 and V^s in the range 
determined above, since they are at most of the order of 0.0001, compared to the leading term 
0.4961. 
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The measured rates may be deduced from the data given by the Particle Data Group |T2 



exp 



2.561 ±0.032 
4.937 ±0.052, 



where, as above, we express all rates in units of 10 MeV; consequently, 

T[K+ ^7f^e+ul,, 



0.519 ±0.016. 



(3.20) 



(3.21) 



exp 



Consequently, the experimental value for R'^^ is 1.023 ± 0.016. 

Up to and including terms of order 0{p^), the theoretical expression for R'^^ is given by 
Eq. (|2.8| ); the terms proportional to m^^^(r — 1)^ and to m^C|' are of order (9(p^). We first 
estimate the 0{p^) correction terms in Eq. ( p.l9| ) by NDA. For the first one, we obtain 







8 r^m^p 


8m3(A^) 


Ml 


NDA 


Ml 


NDA^ M'Ah 



3.2, 



(3.22) 



where p characterizes the size of the >C(o,3) low-energy constants pa, a = 1,2, 3. This implies 
(using = 43.5, from reference 



m^XJr) 



0.22 



i?(M2 - Ml) R 
Similarly, the last term in Eq. ( p.l9| ) is 



0.005. 



(3.23) 



V3R{M^-Ml)\Fl 



r^-1 



^GMO 1/ X 

^- + 3(r2-r) 



r ±2 
r ± 1 



m! 



-0.033 0.006(r2 -r) /r ± 2 
R R 

-0.00075 ± 0.00014 (r2 -r 



r ± 1 
r ±2 

r ± 1 



(3.24) 



which is quite negligible. 

Finally, we turn to the C(p^) corrections to From 



we obtain 



m2^2(r - if 



F3 



- 1 



0.24, 



m i [r - 1) e 

7! 



0.14 5 . 



(3.25) 



(3.26) 
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The NDA estimate for the C2 term gives 



m 



fr- 11 



\_ 2£ 



0.04 ■ 



1 tL 

R ^ Vs. 



(3.27) 



where we note that no assumption has been made here of any suppression of C2 due to Zweig 
rule violation. We have verified that 0{p^) corrections to ei + £2 are negligibly small. Using 
the numerical form of Eq. (|2.5| ), 



o 1 

'=R 



0.533 + 0.502 



1 



.1-- 
IV r2 



(3.28) 



the upshot is that R'^^ may be written (recall Eq. ( p.8| 



R 



+0 



1 + VS £ -0.14 £ ± 



0.22 



(3.29) 



where the term —0.14 e comes from Eq. ( |3.26| ); the indicated uncertainty, coming from other 
higher-order corrections, is dominated by the estimation of the contributions from the -^(0,3) 
terms (see Eq. ( p.23| ) above). Numerically, this implies 



R 



+0 



1 + 



R 



r — 1 / r 
0.85 + 0.80 1 

r + 1 V r2 



± 



0.22 



(3.30) 



and taking, for example, the commonly accepted value [24| R = 43.5 ± 2.2, 



1.020 + 0.005 + 0.018^^^ — ^(1 



r + 1 



r 
r2 



(3.31) 



which is consistent, for any permissible value of r, with the experimental value 1.023 + 0.016 
which we deduced above. In view of this experimental uncertainty, we will leave for a later 
time the careful investigation of the O(p^) effects mentioned above, together with the analysis 
of electromagnetic corrections (the radiative corrections to R^^ in the standard case have been 



investigated in Ref. [0). Note that some (but not all) of the existing experimental data has 
been published with radiative corrections, but often without mention of how these corrections 
have been implemented [Q. More precise knowledge of R~^^ would be useful in constraining 
the relationship displayed in Eq. ( |3.29| ) between the two quark-mass ratios R and r, thereby 
testing the relevance of G^PT. 
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4 Form factor slopes 

Analyses of K^s data frequently assume a linear dependence 



1 + A+,0 



where, as usual, the scalar form factor is defined as 



frit) = f^it) + 



t 



Ml-Ml 



The parameter X^'^ is identical to that of SxPT, 



A 



ML 



+ 



1 



1, 1 Ml, 5 , 1 M%, 



with 



"'f^' + 2 (.-1)3 + 2 (^) - 3 • 



and 



p2 1 ji Y 



12 LI 

^ 327r2 



2 In ^ + In ■ 



M\ 



K 



+ 3 



On the other hand, the GxPT expression for A^'^, 

7l \ 2 



Af^ 

M2 



where 



+ 



- 11 - -(4 - l]^ -m2(r- 1)^3^ 



3847r2Fo2 
1 



5 - 



3 + 



AMI + 1) 



1 



Ml - Ml 



h2{x) 



Ml + Ml 



ho 



M^ 

• TX ■ 

'Ml' 



3847r2Fo2 

(Ml-Ml\ (Ml + MT 
W,+Ml) [mi -Ml 
1 



^GMO 



(r + 3) 



2(M|, + M2) (r-1). 



3 /1 + x 



2 Vl-x 



+ 



3a;(l + x) 
{1-xf 



hix , 
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differs from the S^PT result. This expression for Aq exphcitly displays the dependence on 
m^{r — 1)^A3 and on r, while its dependence on Vus is implicit in the Fx/F^^ terms. However, 
as we saw above in Figure ^ Fk/F.,^ and m^(r — l^A^ are correlated once we know the rate for 
—>■ TT^e^u. Therefore, for any choice of the quark mass ratio r, Eq. ( [4.6[ ) gives a range of 
values for Aq corresponding to a given region in the {Fk/F,,, m^(r — 1)^^43) plane. For example, 
Eq. ( |3.15 ) defines one region in Figure ^, while Eq. (|3.16| ) defines another. We display in Figure 
Hthe corresponding regions in the (Ao,r) plane. 



. 022 



. 021 



. 02 



0.019 



.018 




10 



15 



20 



25 



30 



Figure 4: The dependence of Aq on r. The darker shaded region indicates the range of values 



Eq. ( |3.15| ) for Aq determined by using Vud from superallowed nuclear beta decay, while the 
lighter shaded region shows the corresponding range Eq. ( p.l6| ) determined by using the PDG 
values for Vus- See text for details. 

This analysis was done assuming that the Zweig-violating parameter Zq=0. The whole 
dependence on Zq is contained (through ( = Zq /Aq) in the parameter e(r). As pointed out in 
pT| , the vacuum stability requirement Bq > implies a upper bound on (, which yields 

2 



e(r) -6(r) |c=o < (l-e(r)|^=o) 



(4.8) 



We find that taking ( = or its maximal allowed value makes a difference of less than 0.0004 
in Ao, for any choice of r. 

In SxPT, including one-loop corrections, the result is 

Ao = 0.017 ±0.004, (4.9) 



where the error is an estimate of the uncertainties due to higher-order contributions. The 
experimental situation remains unclear, in view of the inconsistency between some more recent 
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data [0 and the result (Aq = 0.019 ± 0.004) of the high-statistics experiment of Donaldson et 
al p6 . 



5 Summary and discussion 

In the preceding pages, we have studied the Ki^^ form factors at O(p^) precision within the 
generalized framework of chiral perturbation theory. In this connection, three (related) issues 
have been discussed: the extraction of the CKM matrix element Vus from the experimental 
value of the K^^ branching ratios, the determination of the ratio Fx / F^^ of pseudoscalar decay 
constants, and the prediction for the slope Aq of the scalar form factor. 

In the Leutwyler-Roos analysis the expression of (0) is written as 



/f--(0) = l + A + /2 + ---, (5.1) 

where the contributions f\ and /2 are of order mquark and rr?^^^^-^ respectively. The one-loop 
SxPT contribution /i arises from Goldstone boson loops only, a counterterm contribution at 
this order of the standard counting being forbidden by the AdemoUo-Gatto theorem. This 
leads to a parameter free prediction /i = —0.023. Furthermore, /2, which would arise at chiral 



order 0(p ) in S^PT, has been estimated in |TT| to be —0.016 ± 0.008 by using a model for 



pion and kaon wave functions to compute matrix elements in the infinite momentum frame. 
This leads to a value Vus = 0.2196 ± 0.0023. Assuming unitarity of the CKM matrix in a 
three generation standard model, and using the existing estimates of \Vub\, the determination 
Vud = 0.9741±0.0005 from superallowed nuclear /3-decays leads instead to V^s = 0.2261±0.0022. 
Therefore, unitarity of the CKM matrix can only be restored at the expense of having the S^PT 
two- loop correction /2 at least as large as the one-loop contribution /i. This is far from signaling 
a failure of the chiral expansion in the present case, since /i might be anomalously small, being, 
for instance, suppressed, even compared to /2, in the large-iVc limit (another consequence of 
the AdemoUo-Gatto theorem). In G^PT, the corresponding expansion of ff^ (0) reads 

/f'^'(0) = l + /i + /2 + ---, (5.2) 

where /i collects all O(p^) contributions in the generalized chiral counting, and contains 
^(^quark) terms. The difference between /i and fi 



h-h 



8^2 ^ 



(r- 1)^(1 + ^), (5.3) 
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involves a single low-energy constant from ^(2,2)-, ^3, which would appear only at order O(p^) 
in SxPT. The contribution of the £^{2,1) low-energy constant ^ is determined by the value of the 
ratio Fk/ Fj,. The value of Vud from the nuclear /5-decay data can be accommodated by 

(r - l)2m% = -0.124 ±0.022, — = 1.189 ± 0.012 . (5.4) 

F-K 

The corresponding difference is /i — /i = —0.04 ±0.01, and it represents the value which would 
be required in the Leutwyler-Roos analysis for the two-loop contribution /2 (instead of the 
estimate /2 = —0.016 ±0.008 of [|lll]) in order to maintain the unitarity of the CKM matrix. Of 
course, a confirmation, with a comparable accuracy, from other sources (neutron decay, pion 
/9-decay) of the value of Vud obtained from nuclear /3-decays can only be welcome. 

This determination of ^3 and the decrease in the value of the ratio of decay constants, as 
compared to the number F^/ F^^ = 1.22±0.01 |TT[, is compatible with the present experimental 



information concerning the difference in the K^^ and decay rates, and induces only a mild 
modification in the prediction for the slope of the scalar form factor Aq, which, as a function 
of the quark mass ratio m^/m, varies between 0.0018 and 0.0022, well within the range set by 
the high-statistics experiment of Donaldson et al. . The higher values of Aq obtained by 
some of the more recent experiments are therefore difficult to understand at the theoretical 
level, and cannot be ascribed, within G^PT, to the manifestation of a smaller value of the 
bilinear light quark condensate. 

We thus conclude that the nuclear /3-decay determination of Vud = 0.9741±0.0005 need not 
be in contradiction with the present values of the i^es decay rates and with chiral perturbation 
theory. One should then ask how the corresponding increase of \Vus\ by about 2.5 standard 
deviations would manifest itself in various observables. We have already mentioned that the 
present understanding of hyperon semi-leptonic decays is compatible with the suggested update 
of Vus- The effect on the decay rates should be analyzed separately [|23|. Some effect on 
the extraction of the O(p^) low-energy constants Li, L2 and L3 is to be expected a priori, but 
a precise statement requires a closer analysis. Finally, it is worth mentioning the possible effect 
on hadronic spectral functions which are extracted from the decays r —>■ hadrons ± and 
used for a determination of fundamental QCD parameters [^, |3^, ^ . While the non-strange 
(ud) spectral functions should be only barely affected by an increase of ~ 0.25%, the recently 
published strange (us) spectral functions should be reduced by ~ 4.6%. Consequently, we 



would expect no notable influence on the determination of as{M^) [^, whereas the central 
value of the running strange quark mass mg determined recently |]39|, ^ could increase by 
~ 15 — 20%. The issue certainly deserves a more detailed study. 
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Appendices 



A Expansion and renormalization of the effective La- 
grangian 

In this Appendix, we display the structure of the effective action of G^PT up to order 0{p^). 
For a general discussion of the G^PT expansion framework, we refer the reader to the existing 



literature H, 20, 21 



At leading order, the generalized expansion is described by C^'^\ which was first given in 
Ref. i: 

+Ao{{U+xf + ix'-Uf) + Z^iU-'x + X^Uf (A.5) 
+Z^{U^X-X^U)' + Hoix'-X)} ■ 

The notation is as in Refs. [|7|, |^, except for the consistent removal of the factor 2Bq from x, 
the parameter that collects the scalar and pseudoscalar sources, 

X = s + ip = M + -- - , M = diag(m„, rrid, nis) ■ (A.6) 

In G^PT, the next-to-leading-order corrections are of order (9(p^), and still occur at tree level 
only. They are embodied in = C(2,i) + which reads pO|, ^ 



-^^'^ = \F^{i{D^U+D^^U{x-^U + U+x))+i{D^U+D^U){x-^U + U+x) 
+Pi{ix^Uf + (f/+x)') + P2{{X^U + u+x)x^x) 

+P3(x+t/ - U+X){ix^uy - {U+x?) (A.7) 

+P,{{x-'uf + iu+x?){x-'u + U+X) 

+P5{X^X){X^U + U^X) 

+ p,{x^U - u+x)'{x^u + U+x) + P7{X^U + U+x?) ■ 
The tree-level contributions at order O(p^) are contained in Q 

= £(4,0) + -^(2,2) + £(0,4) + -Bo£(0,2) + -So£(2,l) + -So£(0,3) • (^-8) 

^Contributions from the odd intrinsic parity sector are also present in the effective Lagrangian; at order 
©(p^) they are given by the Wess-Zumino term and are the same for SxPT and for GxPT, so that we do not 
display them here. 
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The part without exphcit chiral symmetry breaking, >C{4,o), is described by the same low-energy 
constants Li, L2, Lg and Lio as in S^PT 0: 

£(4,0) = L^{D^U^D>'U){D,U+D'U)+L2{D^U+D^U){D^U+D^U) 

+L-i{D^U+D^'UDM+D''U) - tLgiF^^D^UD^U^ + F^^D^'U^ D^U) (A.9) 

The new term £(2,2) would count as C(p^) in S^PT, and is given by: 
A2,2) = \f^ {A^{D,U+D^Uix^x + U^XX^U)) 

+A,{D,U+D^^U){x+x) 



D^'XX^)) 



+B^{D^U+D^^U{x+Ux^U + U+xU^x)) 
+B2{D^U+xD^U+x + X^D>.Ux^D^U) 
+B^{U+D,xU+D^X + D^x^UD'^x^U) 
+B,{D,U+D^'U){x^Ux^U + U+xU^x) 
+C^{D,Ux^ + xD^U+){D,Ux^ + x^^f/+) 
+C^{D,x^U + f/+Z}^x)(^/.f/+X + X^D^U) 
+C^{D^X^U + f/+D^'x)(^^X+f/ + U+Df^x) 
+C[{D^Ux^ - xD^U+){D^Ux^ - xD^U+) 
+C^{D,x^U - U+D^x){D,U+x - X^D^'U) 
+C^{D^X^U - U+D'^x){D,X^U - U+D'^x) 
+D^D^,U+D^U{x'^U + U+x)){x'^U + U+x) 
+D''{D,U+D^U{x^U - U+x)){x^U - U+x) 
+ H2{D^X^D^X)} ■ 



(A.IO) 



Notice that the number of counterterms (17) involved in £(2,2) agrees with Refs. ||T6|, [17 
However, in both cases, different bases have been used. Finally, the tree-level contributions 
which behave as C'(^quark) ^^e chiral limit are contained in £(0,4), 

\F^{E,{{x^ur+{u+xr) 

+E2{x^xix^Ux^u + u+xu+x)) 

+E3{x^xU^XX^U) 



(0,4) 
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+Ff{x^Ux^U + U+xU+x)' 

+Fi{x^x{x^u + u^x)){x^u + U+x) 

+Ff((x+f/)^ + (t/+x)^)(x+x) 

+F[{x+Ux^U - U+xU^xf (A.ll) 
+Ff ((x+f/)^ - (U^X?) {X^U - U+x) 

+Fi'{x^x{x^u - u+x)){x^u - U+x) 

+Fi'{{x^Uy + {U^xf){x-'U + f/+x)' 

+F^'{X'X){X''U + U+X? 
+FiP{(x^Uf + {U^xf){x'U - U^x? 

+Fr{x^x){x^U-U^xf 

+Fr((x+f/)' - iu^x?){x''U - u^x){x^u + f/+x) 

+ Hs{x^XX^x)+H,{x^xf} ■ 

Contributions from >C(o,4) would only appear at order 0{p^) in S^PT, which, to the best of our 
knowledge, have not been discussed in the literature. 

The O(p^) loop corrections to the processes studied here involve only graphs with one or 
two vertices from C^"^^: 

^(4) _ ^(4) , ^(4) rA 191 

^lloop ~ -^tadpole ^ ^unitarity ^ K-^-^"^) 

The divergent parts of these one-loop graphs have been subtracted at a scale /i in the same 
dimensional renormalization scheme as described in 0]. Accordingly, the low energy constants 
of >C(4_o), j^{2,2), and £(0,4) stand for the renormalized quantities, with an explicit logarithmic 
scale dependence denotes generically any of these renormalized low-energy constants) 

= X(/x') + ■ ln(/i7/i) . (A.13) 

At order (9(p^), the low-energy constants of C^'^^ and C^^^ also need to be renormalized. The 
corresponding counterterms, however, are of order 0{Bq) and 0{Bq), respectively, and they 
are gathered in the three last terms of Eq. ( |A.8|) : in G^PT, renormalization proceeds order by 
order in the expansion in powers of Bq |2y, ^ . Alternatively, one may think of Eqs. ( [A.5|) and 
(|A.7|) as standing for the combinations £*^^^-|-i?Q£(Q 2) and C^^^ + BoCi^2,i)+ BoC'^^^^^-j, respectively, 
with the corresponding low-energy constants representing the renormalized quantities. The full 
list of /3-function coefficients are tabulated below. 
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X 




A, 

yS 
yP 


55oV3 





1052/3 



Table 1: Subtraction scale dependences of the low energy constants of L^'^^ . 



X 


Fi ■ Tx 




3So 


i 


5o 


pi 


ABo (|Ao + Zl + ZP) 


P2 


4fio (1^0 - Zl + ^ZP) 


P3 


ABo (iAo - IZl + IZP) 


pi 


AB, + IZI - \ZP) 


P5 


AB, (Ao + IZl - IZP) 


PQ 


4^0 {t-)ZI + li^cf) 


P7 


4So {%ZI + IZP) 



Table 2: Subtraction scale dependences of the low energy constants of JC'^^\ 



X 




Li 


3 

32 


L2 


3 
16 


L3 





u 


1 
4 




1 

4 


Hi 


1 
8 



Table 3: Subtraction scale dependences of the low energy constants of >C(4,o). 
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A 




A 

Al 


c( ryS '7P\ 


A 

^2 


u 


As 





A, 


2{zl - ZD 


Bi 


(3Ao - IZi - 2Zl) 


B2 


^{zl + ZE) 


B3 





B, 


(A + zl + ZD 


c! 


(Ao - 2ZD 


CI 





cl 





c[ 


(Ao - 2Zi) 







c[ 





D' 


+ 7^(f ) 


DP 


(A + TZD 


H2 






Table 4: Subtraction scale dependences of the low energy constants of C(2,2)- 
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X 

E2 
E3 



f2 
-^0 



r 



X 



-lAl + 8MZ^ + Z^) + f (Zo^ + Z^f 
-IQA,{ZI - ZP) + 2A{Zlf - 24(Z, 



0^ 

7P\2 



26/12 



f (^0" + z^' 



9 42 I 2 4 75 16 4 7P 4 
Z/i(y + 3-^0^0 ^ 9 

9 42 I 2 4 16 4 4 

Z/1q -|- g/lo^o 3"^O^o ~ 9 

-20{z,'r 



7iiz^r 

+ 20(ZoT 



P\2 



29Zo^Zo^ 



Fr 

Ff 



2Al 



^0 



P\2 



—2Aq + ^-AqZq — ^AqZ^ 



{zi + zl 
s{zlf + n{z^f + 2hzlz^ 

29Zn^Z 



n{z^? 



p 

^0 



Fi' 

pSP 

Fr 



44 4 



9- -0^0 

^A^Zl + I 

-^A,ZP + f 
2A,{Zl + Zo^) - 



mzif - 4(z, 



0^^ 



2QZIZI' 



2QZIZP 



A{Zl f - 77(Z, 



A{ZlY + A{Z^f 



2QZlZ^ 
55Z^Zq 



Ha 



8Ao{Z^ 



Zi) + I 



?,h{Zlf + 35(Zo^)2 - 2ZiZ^ 



^^0 9 



17(^0^)2 + n{z^f 



'^ZqZq 



Table 5: Subtraction scale dependences of the low energy constants of >C(o,4). 
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B Masses and Decay Constants 



For the reader's convenience and for later reference, we provide, in this Appendix, the expres- 
sions of the pseudoscalar decay constants and masses at order C(p^) 0. Introducing the notation 
{P = 7r,K,r]) 

Ml , Ml 



where is the renormahsation scale, we obtain the following formulae for the decay constants 
(in the limit m„ = m^). 



with 



with 



p2 _ p2 



1 + 2m4'^ + 2m" 5).%"^ - 4/^^ - 2/i 



2r(2,2) 



K 



(B.2) 



^(2,1) 



v(2,2) 



e + (2 + r)e 



Ai + -^2 + 2^3 + 5i - ^2 

+ i(2 + r')(A4 + 2Bi) + 2(2 + r)D- 



(B.3) 



f2 



^2 
-^0 



1 + 2m5f^^ + 2mHf^^ 



2^- 



(B.4) 



^(2,1) 
'F,K 



s:(2,2) 



-(l + r)e + (2 + r)e 



_(1 + r2)(Ai + A3 + 5i) + -(^2 + 2A3 - 2^2 
+ ^(2 + r=')(A4 + 254) + (1 + r)(2 + r)D^ 



(B.5) 



(B.6) 



^We also take this opportunity to correct a misprint in an earlier published expression |^ of the decay 
constant Fk- 



P2 
-^0 



1 + 2m5^p;^^ + 2m'5'p;;' - Q^ik 
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f(2,l) 1 



)^;^' = (l + 2r)e + (2 + r)e 



= lil + 2r^){A, + ^A2 + 2As + B,-B2) 

+1(2 + r^)(A^ + 2S4) + ^(1 + 2r)(2 + r)D^ 
-l{r-mC[ + Ci) 
For the masses, we obtain: 

F^M^ = F2|2mBo + 4m% + 4m2(2 + r)Zo^ 

+2m35SJ + 2m^5K + Am'A.M'^ 

- /X. [3 + 5 e(r)] 
-2/.^ [M^ + (l + r)M2e(r)] 

4/^''K + (l + 20M^»-2^ 



+ r)mBo + (1 + r)2m2>lo + 2(1 + r)(2 + r)m2Zo^ 

+(1 + r)m-^5£'i + (1 + r)m^(5£tJ + (1 + r)=^m'^M|: 



5 Mi + (r + 1) (1 + 2r) e(r) + 



r + 1 



+ 2r)mBo + ^(1 + 2r2)m2>lo + ^(1 + 2r)(2 + r)m'Zo^ + ^(1 



+2m3<;J + 2m^5£;y + ^(1 + 2r')rn'A,M'^ - ^(1 - r)2m2crM, 
-^^.[M^^ + i2r + l)M^^eir)-2^' 



1 



- 3M^ + (r + l)(2r + 1) Ml e(r) + 2 ^— ^ Agmo 

16Mi - + 3 (2r + 1)^ Ml e(r) + 4 A^mo 

r — i 
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with 



<^£? = ^Pi + ^P2 + (10 + 4r + r>4+^(2 + r>5 + 6(2 + r)V 



x(o,3) 3 



(B.ll) 

'M,K = -(l + r + r>i + i(l-r + r')p2 + 3(2 + 2r + r>4 + ^(2 + r2)p5 + 6(2 + r)V 

(B.12) 

= ^(l + 2r>i + i(l + 2r>, + ^(l + r)(l-r)V3 + 

^(10 + 8r + ITr^ + 10r3)p4 + ^(1 + 2r)(2 + r^jps + 

^(1 - rfi2 + r)pe + 2(1 + 2r)(2 + r) V , (B.13) 
while the C^(p^) contributions are 

= 8^1 + 2^2 + 4(1 + r) (2 + r^)^^ + (20 + 9r + r-^)Ff 

+ (4 + r + r^)Fi + 2(2 + r^)^/ + 4(2 + r)(6 + 2r + r2)F5^^ (B.14) 
+2(2 + r)(2 + r')F6^^ , 

= 2(l + r)(l + r^)£;i + (l + r=^)£;2 + ^(l + r)(l-rr£;3 
+4(1 - r) (2 + r2)Ff + (8 + 9r + 9r^ + 4:r^)Fi 

+ (4 - r + + 2r^)Fi + (1 + r)(2 + r^)Ff (B.15) 
+4(2 + r)(4 + 3r + 2r2)F5^^ + 2(2 + r)(2 + r2)F6^^ , 

C = ^(l + 2r^)£;, + |(l + 2r^)£;2 



+^(1 + 2r2)(2 + r2)Ff + J(20 + 13r + 37r^ + 20r^)F, 

+^(4 + 5r + Sr^ + 4:r^)F^ + ^(1 + 2r2)(2 + r^)F^ 

+ — (1 + rf{l - rfF[ + 4(1 - rf{l + r + r^)Fi 
3 



5 



(1 - rf{l + r + r^)F3^ + 4(2 + r)(2 + 2r + 3r^ + 2r^)F5^^ 



+ |(2 + r)(2 + r')(l + 2r)F6^^ + ^(1 - rf{2 + r'jF^^^ (B.16) 
+^(1 - r)2(2 + r2)F6^^ + ^(1 - 0^(1 + r)(2 + ^F^^^ . 
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